SSE: a nucleotide and amino acid sequence analysis platform

<p>Abstract</p> <p>Background</p> <p>There is an increasing need to develop bioinformatic tools to organise and analyse the rapidly growing amount of nucleotide and amino acid sequence data in organisms ranging from viruses to eukaryotes.</p> <p>Finding</p> <p>A simple sequence editor (SSE) was developed to create an integrated environment where sequences can be aligned, annotated, classified and directly analysed by a number of built-in bioinformatic programs. SSE incorporates a sequence editor for the creation of sequence alignments, a process assisted by integrated CLUSTAL/MUSCLE alignment programs and automated removal of indels. Sequences can be fully annotated and classified into groups and annotated of sequences and sequence groups and access to analytical programs that analyse diversity, recombination and RNA secondary structure. Methods for analysing sequence diversity include measures of divergence and evolutionary distances, identity plots to detect regions of nucleotide or amino acid homology, reconstruction of sequence changes, mono-, di- and higher order nucleotide compositional biases and codon usage.</p> <p>Association Index calculations, GroupScans, Bootscanning and TreeOrder scans perform phylogenetic analyses that reconcile group membership with tree branching orders and provide powerful methods for examining segregation of alleles and detection of recombination events. Phylogeny changes across alignments and scoring of branching order differences between trees using the Robinson-Fould algorithm allow effective visualisation of the sites of recombination events.</p> <p>RNA secondary and tertiary structures play important roles in gene expression and RNA virus replication. For the latter, persistence of infection is additionally associated with pervasive RNA secondary structure throughout viral genomic RNA that modulates interactions with innate cell defences. SSE provides several programs to scan alignments for RNA secondary structure through folding energy thermodynamic calculations and phylogenetic methods (detection of co-variant changes, and structure conservation between divergent sequences). These analyses complement methods based on detection of sequence constraints, such as suppression of synonymous site variability.</p> <p>For each program, results can be plotted in real time during analysis through an integrated graphics package, providing publication quality graphs. Results can be also directed to tabulated datafiles for import into spreadsheet or database programs for further analysis.</p> <p>Conclusions</p> <p>SSE combines sequence editor functions with analytical tools in a comprehensive and user-friendly package that assists considerably in bioinformatic and evolution research.</p

LAVA: An Open-Source Approach To Designing LAMP (Loop-Mediated Isothermal Amplification) DNA Signatures

<p>Abstract</p> <p>Background</p> <p>We developed an extendable open-source Loop-mediated isothermal AMPlification (LAMP) signature design program called LAVA (LAMP Assay Versatile Analysis). LAVA was created in response to limitations of existing LAMP signature programs.</p> <p>Results</p> <p>LAVA identifies combinations of six primer regions for basic LAMP signatures, or combinations of eight primer regions for LAMP signatures with loop primers, which can be used as LAMP signatures. The identified primers are conserved among target organism sequences. Primer combinations are optimized based on lengths, melting temperatures, and spacing among primer sites. We compare LAMP signature candidates for <it>Staphylococcus aureus </it>created both by LAVA and by PrimerExplorer. We also include signatures from a sample run targeting all strains of <it>Mycobacterium tuberculosis</it>.</p> <p>Conclusions</p> <p>We have designed and demonstrated new software for identifying signature candidates appropriate for LAMP assays. The software is available for download at <url>http://lava-dna.googlecode.com/</url>.</p

Correlation between nucleotide composition and folding energy of coding sequences with special attention to wobble bases

Background: The secondary structure and complexity of mRNA influences its accessibility to regulatory molecules (proteins, micro-RNAs), its stability and its level of expression. The mobile elements of the RNA sequence, the wobble bases, are expected to regulate the formation of structures encompassing coding sequences. Results: The sequence/folding energy (FE) relationship was studied by statistical, bioinformatic methods in 90 CDS containing 26,370 codons. I found that the FE (dG) associated with coding sequences is significant and negative (407 kcal/1000 bases, mean +/- S.E.M.) indicating that these sequences are able to form structures. However, the FE has only a small free component, less than 10% of the total. The contribution of the 1st and 3rd codon bases to the FE is larger than the contribution of the 2nd (central) bases. It is possible to achieve a ~ 4-fold change in FE by altering the wobble bases in synonymous codons. The sequence/FE relationship can be described with a simple algorithm, and the total FE can be predicted solely from the sequence composition of the nucleic acid. The contributions of different synonymous codons to the FE are additive and one codon cannot replace another. The accumulated contributions of synonymous codons of an amino acid to the total folding energy of an mRNA is strongly correlated to the relative amount of that amino acid in the translated protein. Conclusion: Synonymous codons are not interchangable with regard to their role in determining the mRNA FE and the relative amounts of amino acids in the translated protein, even if they are indistinguishable in respect of amino acid coding.Comment: 14 pages including 6 figures and 1 tabl

Unique and conserved MicroRNAs in wheat chromosome 5D revealed by next-generation sequencing

MicroRNAs are a class of short, non-coding, single-stranded RNAs that act as post-transcriptional regulators in gene expression. miRNA analysis of Triticum aestivum chromosome 5D was performed on 454 GS FLX Titanium sequences of flow sorted chromosome 5D with a total of 3,208,630 good quality reads representing 1.34x and 1.61x coverage of the short (5DS) and long (5DL) arms of the chromosome respectively. In silico and structural analyses revealed a total of 55 miRNAs; 48 and 42 miRNAs were found to be present on 5DL and 5DS respectively, of which 35 were common to both chromosome arms, while 13 miRNAs were specific to 5DL and 7 miRNAs were specific to 5DS. In total, 14 of the predicted miRNAs were identified in wheat for the first time. Representation (the copy number of each miRNA) was also found to be higher in 5DL (1,949) compared to 5DS (1,191). Targets were predicted for each miRNA, while expression analysis gave evidence of expression for 6 out of 55 miRNAs. Occurrences of the same miRNAs were also found in Brachypodium distachyon and Oryza sativa genome sequences to identify syntenic miRNA coding sequences. Based on this analysis, two other miRNAs: miR1133 and miR167 were detected in B. distachyon syntenic region of wheat 5DS. Five of the predicted miRNA coding regions (miR6220, miR5070, miR169, miR5085, miR2118) were experimentally verified to be located to the 5D chromosome and three of them : miR2118, miR169 and miR5085, were shown to be 5D specific. Furthermore miR2118 was shown to be expressed in Chinese Spring adult leaves. miRNA genes identified in this study will expand our understanding of gene regulation in bread wheat

Dynamic probe selection for studying microbial transcriptome with high-density genomic tiling microarrays

<p>Abstract</p> <p>Background</p> <p>Current commercial high-density oligonucleotide microarrays can hold millions of probe spots on a single microscopic glass slide and are ideal for studying the transcriptome of microbial genomes using a tiling probe design. This paper describes a comprehensive computational pipeline implemented specifically for designing tiling probe sets to study microbial transcriptome profiles.</p> <p>Results</p> <p>The pipeline identifies every possible probe sequence from both forward and reverse-complement strands of all DNA sequences in the target genome including circular or linear chromosomes and plasmids. Final probe sequence lengths are adjusted based on the maximal oligonucleotide synthesis cycles and best isothermality allowed. Optimal probes are then selected in two stages - sequential and gap-filling. In the sequential stage, probes are selected from sequence windows tiled alongside the genome. In the gap-filling stage, additional probes are selected from the largest gaps between adjacent probes that have already been selected, until a predefined number of probes is reached. Selection of the highest quality probe within each window and gap is based on five criteria: sequence uniqueness, probe self-annealing, melting temperature, oligonucleotide length, and probe position.</p> <p>Conclusions</p> <p>The probe selection pipeline evaluates global and local probe sequence properties and selects a set of probes dynamically and evenly distributed along the target genome. Unique to other similar methods, an exact number of non-redundant probes can be designed to utilize all the available probe spots on any chosen microarray platform. The pipeline can be applied to microbial genomes when designing high-density tiling arrays for comparative genomics, ChIP chip, gene expression and comprehensive transcriptome studies.</p

Selective Light-Triggered Release of DNA from Gold Nanorods Switches Blood Clotting On and Off

Blood clotting is a precise cascade engineered to form a clot with temporal and spatial control. Current control of blood clotting is achieved predominantly by anticoagulants and thus inherently one-sided. Here we use a pair of nanorods (NRs) to provide a two-way switch for the blood clotting cascade by utilizing their ability to selectively release species on their surface under two different laser excitations. We selectively trigger release of a thrombin binding aptamer from one nanorod, inhibiting blood clotting and resulting in increased clotting time. We then release the complementary DNA as an antidote from the other NR, reversing the effect of the aptamer and restoring blood clotting. Thus, the nanorod pair acts as an on/off switch. One challenge for nanobiotechnology is the bio-nano interface, where coronas of weakly adsorbed proteins can obscure biomolecular function. We exploit these adsorbed proteins to increase aptamer and antidote loading on the nanorods.National Science Foundation (U.S.) (Grant DMR #0906838

The concept of RNA-assisted protein folding: the role of tRNA

We suggest that tRNA actively participates in the transfer of 3D information from mRNA to peptides - in addition to its well-known, "classical" role of translating the 3-letter RNA codes into the one letter protein code. The tRNA molecule displays a series of thermodynamically favored configurations during translation, a movement which places the codon and coded amino acids in proximity to each other and make physical contact between some amino acids and their codons possible. This specific codon-amino acid interaction of some selected amino acids is necessary for the transfer of spatial information from mRNA to coded proteins, and is known as RNA-assisted protein folding

DNA Suspension Arrays: Silencing Discrete Artifacts for High-Sensitivity Applications

Detection of low frequency single nucleotide polymorphisms (SNPs) has important implications in early screening for tumorgenesis, genetic disorders and pathogen drug resistance. Nucleic acid arrays are a powerful tool for genome-scale SNP analysis, but detection of low-frequency SNPs in a mixed population on an array is problematic. We demonstrate a model assay for HIV-1 drug resistance mutations, wherein ligase discrimination products are collected on a suspension array. In developing this system, we discovered that signal from multiple polymorphisms was obscured by two discrete hybridization artifacts. Specifically: 1) tethering of unligated probes on the template DNA elicited false signal and 2) unpredictable probe secondary structures impaired probe capture and suppressed legitimate signal from the array. Two sets of oligonucleotides were used to disrupt these structures; one to displace unligated reporter labels from the bead-bound species and another to occupy sequences which interfered with array hybridization. This artifact silencing system resulted in a mean 21-fold increased sensitivity for 29 minority variants of 17 codons in our model assay for mutations most commonly associated with HIV-1 drug resistance. Furthermore, since the artifacts we characterized are not unique to our system, their specific inhibition might improve the quality of data from solid-state microarrays as well as from the growing number of multiple analyte suspension arrays relying on sequence-specific nucleic acid target capture

Assessment and validation of a suite of reverse transcription-quantitative PCR reference genes for analyses of density-dependent behavioural plasticity in the Australian plague locust

<p>Abstract</p> <p>Background</p> <p>The Australian plague locust, <it>Chortoicetes terminifera</it>, is among the most promising species to unravel the suites of genes underling the density-dependent shift from shy and cryptic solitarious behaviour to the highly active and aggregating gregarious behaviour that is characteristic of locusts. This is because it lacks many of the major phenotypic changes in colour and morphology that accompany phase change in other locust species. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) is the most sensitive method available for determining changes in gene expression. However, to accurately monitor the expression of target genes, it is essential to select an appropriate normalization strategy to control for non-specific variation between samples. Here we identify eight potential reference genes and examine their expression stability at different rearing density treatments in neural tissue of the Australian plague locust.</p> <p>Results</p> <p>Taking advantage of the new orthologous DNA sequences available in locusts, we developed primers for genes encoding 18SrRNA, ribosomal protein L32 (RpL32), armadillo (Arm), actin 5C (Actin), succinate dehydrogenase (SDHa), glyceraldehyde-3P-dehydrogenase (GAPDH), elongation factor 1 alpha (EF1a) and annexin IX (AnnIX). The relative transcription levels of these eight genes were then analyzed in three treatment groups differing in rearing density (isolated, short- and long-term crowded), each made up of five pools of four neural tissue samples from 5^thinstar nymphs. SDHa and GAPDH, which are both involved in metabolic pathways, were identified as the least stable in expression levels, challenging their usefulness in normalization. Based on calculations performed with the geNorm and NormFinder programs, the best combination of two genes for normalization of gene expression data following crowding in the Australian plague locust was EF1a and Arm. We applied their use to studying a target gene that encodes a Ca²⁺binding glycoprotein, <it>SPARC</it>, which was previously found to be up-regulated in brains of gregarious desert locusts, <it>Schistocerca gregaria</it>. Interestingly, expression of this gene did not vary with rearing density in the same way in brains of the two locust species. Unlike <it>S. gregaria</it>, there was no effect of any crowding treatment in the Australian plague locust.</p> <p>Conclusion</p> <p>Arm and EF1a is the most stably expressed combination of two reference genes of the eight examined for reliable normalization of RT-qPCR assays studying density-dependent behavioural change in the Australian plague locust. Such normalization allowed us to show that <it>C. terminifera </it>crowding did not change the neuronal expression of the <it>SPARC </it>gene, a gregarious phase-specific gene identified in brains of the desert locust, <it>S. gregaria</it>. Such comparative results on density-dependent gene regulation provide insights into the evolution of gregarious behaviour and mass migration of locusts. The eight identified genes we evaluated are also candidates as normalization genes for use in experiments involving other Oedipodinae species, but the rank order of gene stability must necessarily be determined on a case-by-case basis.</p